Anticipatory synchronization with variable time delay and reset 
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A method to synchronize two chaotic systems with anticipation or lag, coupled in the drive 
response mode is proposed. The coupling involves variable delay with three time scales. The method 
has the advantage that synchronization is realized with intermittent information about the driving 
system at intervals fixed by a reset time. The stability of the synchronization manifold is analyzed 
with the resulting discrete error dynamics. The numerical calculations in standard systems like the 
Rossler and Lorenz systems are used to demonstrate the method and the results of the analysis. 
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I. INTRODUCTION 



The synchronization of uni-directionally coupled 
chaotic systems have been studied reasonably well in the 
past few years [l], [E [I[ ■ The synchronization state in such 
cases can be phase, lag, generalized or complete depend- 
ing upon the strength of the coupling [J, |5|. Recently 
synchronization of systems via coupling with a time de- 
lay, which presumably takes care of the finite propaga- 
tion times, switching speeds and memory effects have 
been reported [E 0, IE IE] ■ Such studies relate to a vari- 
ety of diverse phenomena like chirping of crickets, neu- 
ral networks, automatic steering and control and cou- 
pled phase locked lasers [IE El, EE] • When the coupling 
is not isochronous with the system dynamics, it is possi- 
ble to realize retarded (delay) , complete and anticipatory 
synchronization of chaotic systems. Moreover, synchro- 
nization in such cases reveal many novel phenomena like 
parametric resonance 13j , multi-stable phase clustering 
Ul EE EE E3, amplitude death etc QJ, EE S3- An 
interesting aspect of such delay induced synchronization, 
that has attracted lot of attention, is that the driven sys - 
tcm can anticipate the dynamics of the driver [U [2J] • 
The maximum possible anticipation time is reported to 
be enhanced considerably by using an array or ring of 
such systems [IE HE HE)]. Experimental verification of 
anticipatory and retarded synchronization is reported in 
electronic circuits as well as semiconductor lasers with 
delayed optoelectronic feed back [IE l2Sj |- In all these 
studies, the delay time in the coupling, once chosen, re- 
mains constant as the system evolves. 

The synchronization of chaotic systems, in general, has 
attracted great attention due itspotential application in 
secure communication [29l l30l l3lj |. However the use of 
low dimensional systems in this context is found to be 
insecure due to the ease of reconstruction from the trans- 



mitted signal [32, |33j ■ Therefore recently, chaos synchro- 
nization in high dimensional systems, especially systems 
with an inherent time delay, has been proposed as a bet- 
ter alternative 0, HE HE M, EH ■ 

In this paper we propose a method of de- 
lay/anticipatory synchronization with coupling involving 
variable time delay. Here, the synchronization can be re- 
alized with limited information about the driver via oc- 
casional contacts or feedbacks at specific intervals. This 
makes the method highly cost effective and can be ap- 
plied to cases where the signal transmission from driver is 
slow or intermittent. This is achieved by using a variable 
delay in the coupling that is reset at definite intervals. 
The dynamics then evolves under three additional time 
scales, the delay t±, the anticipatory time r 2 and the re- 
set time r. Unlike the case of fixed delay, the resetting 
mechanism makes the error dynamics discrete and it is 
possible to carry out an approximate analytic analysis. 
The analysis gives the maximum r 2 for a given r. This 
also fixes the regions of stability in the parameter plane 
of coupling and delay. The method is demonstrated for 
standard systems like Rossler and Lorenz. 



II. SYNCHRONIZATION WITH VARYING 
DELAY AND RESET 

A. Model system 

Consider a dynamical system x of dimension n that 
drives an identical system y. We choose a simple cou- 
pling term of the linear difference type but with the drive 
variable delayed by t\ and the driven variable delayed by 
t%. Thus, the dynamics is given by 

x = f(x) (la) 

OO 

y = f(y) + e^2r(x tl -y t2 )x(mT,( m +i)T) (it>) 
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where x tl = x(t — t{), y t2 = y(t — i 2 ), r is the reset- 
ting time and X(t',t") is an indicator function such that 
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X(t',t") = 1 for t' < t < t" and zero otherwise. Here, 
r = \T\, T2, . . • , r„] T is a constant vector specifying the 
coupling between the components of x and y. In nu- 
merical simuations we take only one coponent of T to be 
nonzero. Both the delays t\ and £2 depend on time and 
we choose this dependence as 

U = n + t — mr, i = 1,2. 

Thus, t — ti = ttit — 7*. As the two systems evolve, t\ and 
ti also evolve with the same time scale and the coupling 
term uses the same value of both variables x^ and y t2 
during each resetting time interval r, i.e. the coupling 
term is constant for the time interval r. In each time 
interval r, the initial values of the delays t\ and ti are 
T\ and r 2 respectively. The delays increase linearly with 
time upto values ti+t and T2 + r and then they are reset 
for the next interval. As a consequence, the coupling 
requires the variable of the drive system only at discrete 
time intervals of r. We also note that t\ — ti = n — Ti 
for all t. 



B. Synchronization Manifold 

Synchronization manifold for the coupled systems (TTJ) 
is defined by y(t — r 2 ) = x(t — n) or y(t) = x(t — T\ +J2) 
Thus, we can get all the following three possibilities [26j |. 
(1) If T\ — T 2 > 0, we can get delay or lag synchronization 
with T\ — r 2 as the lag time. (2) If t\ — r 2 < 0, we 
can get anticipatory synchronization with T2 — t\ as the 
anticipation time. (3) If T\ — t 2 — 0, we can get equal 
time synchronization. 

As an illustration of this, we take the standard Rossler 
oscillator in the chaotic state as the driver described by 
the equations 

x\ = -x 2 - x 3 

x\ = x\ + ax2 (2) 
X3 = b + 2:3(2:1 - c). 

This is coupled to an identical system through the cou- 
pling scheme given in Eq. (fTJ). Only x\ and y\ are cou- 
pled, i.e. T = [1,0, 0] T . Taking the parameter val- 
ues a = 0.15,6 = 0.2 and c = 10.0, both the systems 
are evolved from random initial conditions using Runge 
Kutta algorithm with a time step 0.01 for 2000 units of 
time. With r = 0.10 and the coupling strength e = 0.4, 
the resulting time series obtained for Ti=0.84 and T2 
= 0.02 is plotted in Fig.l.a. Here the response system 
y(t) (dashed line) lags behind the driver x(t) (solid line) 
by T2 — Ti. Fig.l.b shows the same for t\= 0.02 and 
T2=0.84 where y(t) anticipates x(t) with the same time 
shift. The degree of synchronization with the correspond- 
ing time shift can be quantified using the similarity func- 
tion defined as 

S 2 (T) = < [yi(t)- Xl (t + T)] 2 > 

V< xW) >< VW) > 




0.8 0.82 0.84 0.8 0.82 0.84 



T T 

FIG. 1: The simulated time series of two Rossler systems 
coupled through the scheme in Eq. |l]). In (a) the case of delay 
synchronization is shown with a delay of 0.82 units between 
the xi(t) (solid line) and j/i (t) (dashed line), (b) is a case of 
anticipatory synchronization when yi(t) anticipates Xi(t) by 
the same units. The similarity function S 2 (T) corresponding 
to both these cases are shown in (c) and (d) respectively. 

Figs.l.c and l.d show S 2 (T) computed for different values 
of T. The minimum occurs at 0.82, i.e. T — \t\ — r 2 |, 
indicating synchronization with delay or anticipation of 
0.82 time units. 

It should be noted that the delay time t\ is not of much 
significance in the error dynamics, since the time scale of 
the drive system can be linearly shifted by t± . This point 
will become clear when we do the stability analysis in the 
next section. 



III. LINEAR STABILITY ANALYSIS 

The dynamics of the system in Eq. (jTJ) involves three 
time scales in addition to its inherent scale. Define the 
transverse system by the variable A = y — x Tl - T2 . Its 
dynamics in linear approximation can be derived from 
Eq. $TJ) as 

OO 

A = / (x Tl _ T2 )A - £ ^ X(mr,(m+l)T)Am (4) 
m=0 

where A m = A(i— £2) = A(mr— T2) and we take coupling 
in all components of x and y, i.e. T = [1, 1 . . . , 1] T . Thus, 
A m is a constant in each time interval mr < t < (m + 
1)t. We note that t\ enters only through the Jacobian 
term / and can be eliminated by shifting the time scale 
of the drive system linearly and redefining T2 suitably. 
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Hence, as noted in the previous section, n is not very 
significant for the stability analysis. The fixed point A = 
corresponds to the lag/anticipatory synchronized state. 

In general, it is not possible to solve Eq. (fj|. However, 
we can approximate the equation by replacing Jacobian 
/ by some effective time average Lyapunov exponent A 
(only the real part is required). 



with the solutions 



A = AA - £ ^2 X(mr,(m+l)r)A m 
rn— 



(5) 



In the following analysis we assume A to be positive. The 
results can be easily extended to A < (see Appendix [C]) . 

^From the numerical analysis presented in the next 
section, it appears that the approximation of replacing 
/ by an effective A is reasonable for small values of T2. 
We need a larger value of A for large r 2 . 

In the interval tot < t < (to + l)r, the solution of 
Eq. © is 



A = oA m + C m e 



(0) 



where a = e/A is the normalized dimensionless coupling 
constant, and C m is an integration constant. 



A. < r 2 < r 

Let us first consider the case < t 2 < r. For t — (to + 
1)t — t 2 , A = A m+1 . Thus, eliminating the integration 
constant, Eq. © gives 

A = oA m + (A m+1 - aA m )e x ^ m+1 ^+^ (7) 

For (m — l)r < t < tot we have 

A = aA m _! + (A m - aA m _ 1 )e A ( t - raT+ ") (8) 

Matching the solutions (O and ((SJ at t — tot, and sim- 
plifying, we get the following recursion relation 

A m+1 = a(l - e x ^- T ^ + -e Xr )A m 

a 

-ae Xr (l - e- AT3 )A m _! (9a) 
= aA m -6A m _! (9b) 



where 



b = ae XT {l-e- XT2 ). 



(10a) 
(10b) 

We can write Eq (|9b[) as a 2-d map in matrix form as, 
A ™+i \ _ fa -b \ f A m \ , , 

A m J ~ \ 1 ) \ A m _ x ; [n > 

The eigenvalue equation for the Jacobian matrix is 

H 2 -an + b = (12) 



M± = \( a± Va 2 - 46) 



(13) 



The synchronized state, A = 0, is stable if both the so- 
lutions satisfy \fj,±\ < 1. The detailed analysis of the 
stability conditions is given in Appendix [SJ 

Fig. [2] shows the stability region in t 2 /t — a plane. 
The lower limit of stability is always ai — 1. For smaller 
values of t 2 (t 2 < t 2j ,), the upper limit of stability is 
given by (Eq. (|A6[0 



(14) 



2e A ( r " T 2) - e Ar - 1 



while for larger values of r 2 (r 2p < t 2 < r) it is given by 
(Eq. {AS])) 



,-At 



1 - e- Ar 2 



(15) 



The maximum value of a p is given by the intersection of 
the two curves ([H]) and (fT5|) . 



3 + e AT 

The corresponding r 2p value is given by 

a p (a p + 3) 



T 2 p 



(a„ + l) 2 



(16) 



(17) 



We also obtain r 2ma:I ., the maximum allowed value of 
t 2 for the stability of the synchronized state (Eq. (|A9j) ). 
A general expression for r 2 „ laa; is obtained in the next 
subsection (Eq. [[29"])), 



B. T2 > T 

Let t 2 = fcr + r 2 , k = 0, 1, . . . where r 2 < r. Consider 
the solution ([6]) in the interval mr < t < (m + l)r. Then 
for t = (to + 1)t — r 2 = (m + fc + l)r — t 2 , we get 

A m+k+1 = aA m + C m e x ^ +1 ^- x < (18) 
Hence Eq. ^ becomes 

A = aA m + (A m+fc+1 - aA m )e x ^ m+1 ^ + ^ (19) 
For (to — l)r < t < mr we have 

A = aA ro _! + (A m+fc - aA m _ 1 )e A ( t - mT +^) (20) 

Equating the solutions l[20|) and p9|) for t = tot, and 
simplifying we get the following recursion relation 

A m+k+1 = e Xr A m+k - a(e x( - T - T '^ - l)A m 

-ae XT (l - e- A ^)A m _! (21) 
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FIG. 2: This figure shows the stability region of the syn- 
chronized state (A = solution of Eq. (O) in the T2/V — a 
plane. The solid line is for Ar = 0.25 and the dashed line 
is for Ar = 0.5. The lower limit of stability is a; = 1 (dot- 
ted line). For smaller values of r 2 < r 2p , the upper limit 
of stability is given by (Eq. (|14[l) while for larger values of 
t~2 (r 2p < T2 < r) it is given by (Eq. (fl5jl 'l. The peak val- 
ues are (0.0530 . . . /0.25 . . . = 0.212, 15.083 . . .) for Ar = 0.25 
and (0.088... /0.5 = 0.176 7.169 .. .)for Ar = 0.5 (see 
Eqs. ((T6j and (JT7J). For r < r 2 < 2r the upper limit of 
stability is given by Eq. (|27p .The maximum value of Ar 2 is 
0.74 for Ar = 0.5. 



The lower stability limit is oti = 1. The upper stability 
limit can be obtained by noticing that when the magni- 
tude of the imaginary [i becomes one, the two imaginary 
eigenvalues can be written as /i = e ±l6 and the above 



equation has a factor /i 2 
condition 



2cos(9)fi + 1. This gives the 



cb o + b 1 -l = (25) 
Using this we get a quadratic equation for a. 



a 2 a 



a\a — 1 = 



(26) 



where a\ = e 2AT (l - e~ Ar 2) + e A(r ~ T 2) _ i and a 2 = 

e 2Ar (l — e~ AT s>) 2 . One solution of this equation gives the 
upper stability limit for a. 



1 



(-ai + Jal+4a 2 ) 



(27) 



This upper stability limit is shown in Fig. [3] for r < r2 < 
2r. 

We can also obtain T 2max , the maximum value of r 2 for 
which synchronization is possible. This happens when 
there is always an eigenvalue with magnitude greater 
than one, i.e. when ai — a u — 1. By putting a u = ai — 1 
in Eq. (|26|) . we get ai + a 2 = 1. However a better condi- 
tion is obtained if we note that for a u = a/ = 1, Eq. (12"4")) 
has two degenerate solutions /z = 1, i.e. Eq. (|2"4"|) has a 
factor /i 2 — 2/i + 1. This gives the conditions 60 = 2 — c 
and 1 — bi = 2(2 — c). First condition gives 



This gives a map of dimension k + 2. In matrix form, the 



map can be expressed as 



A. 




1 
1 



h 60 \ 



/ A m+ fe + i \ 



A m _(_fc 
A m _|_fc_i 



\00 ... 1 0/\A 



where c = e Ar , 

-Ar.' 



61 = a(e x( - T - T ^ - 1) and b 



e 2 ). The eigenvalue equation is 



f i k+2 - c f i k+1 + b lf x + b Q = 



(22) 
ae AT (l - 

(23) 

For k = 0, the map of Eq. (|2"2"|) reduces to the 2d-map of 
Eq. (fTTj) . In general the behavior of the largest magni- 
tude [i is as shown in Fig. [TJa, i.e. the stability range is 
from ai = 1 till the complex fi has magnitude one. 

1. k = 1, i.e. r < r 2 < 2r 



Ar 2r . 



Ar + Ar 2ma3 . 
Ar-ln2-ln(l 



2. General k 



-Xt\ 



(28) 



For a general getting explicit solutions for a u is not 
easy. But it is possible to get an expression for T 2max . 
We use the condition that for a u = a.\ = 1, Eq. (f2"3")l has 
two degenerate solutions /1 = 1 . This gives the condition 
60 = 1 + k — kc. Simplifying, we get 



k\r + \T 2max 

k\r - ln(fe + 1) - ln(l - e~ AT ) (29) 



Ar 2r , 



For k — 0, this equation reduces to Eq. (|A9|) and for 
k = 1 it reduces to Eq. (|28| . Fig. [3] shows the plot of 
X.T 2ma x as a function of r for different k values. For each 
k the plot is for the range (t^, Tk+i) where r^ is defined by 
krk = T 2max and from Eq. (|29p we get rfe = ln((/c + l)/fc). 



For k = 1, we have a 3d-map. The eigenvalue equation 
becomes 



/i 3 - c/i 2 + bin + b = 



(24) 



r 2 



In this case it is not possible to obtain the stability 
range for the synchronized state in terms of a. However, 
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0.2 0.4 0.6 0.8 1 
At 



FIG. 3: The figure plots the maximum T2max as a func- 
tion of At (solid line). Here, T2max, k = 0, 1, 2, . . . is given 
by Eq. (|29(l ) and for each k the range of r is (rk,Tk+i). 
The dashed curve passes through the values T2max = kr^ = 
fcln((fc+l)/fc) (Eq. (|30|l ). The inset shows the same plot with 
Ar range (0,5). 

it is possible to obtain an explicit expression for Timax as 
(The detailed calculations are given in Appendix |B1) 

^T 2max = Ar/(exp(Ar) - 1) (30) 

The dashed line in Fig. ([3]) corresponds to Eq. (|30|) . It 

gives the correct Timax only for t 2 = nr. 

IV. NUMERICAL ANALYSIS 

We choose two standard systems, Rossler and Lorenz, 
to confirm the main results obtained in the previous sec- 
tion. 

In the case of Rossler system given in Eq. ([2]), two 
identical systems are coupled in the drive response mode 
via the coupling scheme explained in Eq. JTJ). Start- 
ing from random initial conditions and choosing the sys- 
tem parameters in the chaotic region (a = 0.15,6 = 0.2 
and c = 10.0), they are evolved for 200000 units with 
a time step of 0.01. The correlation coe fficient C =< 
yi(t)xi{t + t 2 ) > /yj< x\{t) >< yj(t) > between x x (t) 
and 2/1 (t) shifted by the effective r 2 — |t 2 — t\\ (hereafter 
referred to as t 2 itself) is calculated using the last 5000 
values. The region of stability of the synchronized state 
is isolated as the region where C = 0.99 and boundaries 
of stability fixed when C goes below this value. 

Taking r = 0.5, r 2 is varied from to 1.0 units in steps 
of 0.01. For each value of r 2 , the coupling strength e 
is increased in steps of 0.005. The appropriately shifted 
correlation coefficient is calculated and using the criterion 
mentioned above the lower and upper limits of stability 
are found out. The results are plotted in the parameter 




0.2 0.4 0.6 0.8 1 

FIG. 4: The limits of stability of the synchronized state of two 
chaotic Rossler systems in the parameter plane T2 — e. The 
solid line is the limits obtained from the stability analysis for 
A = 0.65. The agreement is good for values of e > 3.0. 

plane t 2 — e in Fig.|4j The overall behavior agrees with the 
theoretical analysis carried out in the previous section. 
The upper limits obtained by the stability analysis given 
in Eq. (fli]). Eq. (JT5J) and Eq. ((27]) for the different relative 
ranges of r 2 are calculated for a typical value of A = 0.65 
and shown as solid line. For values of r 2 < 0.3 or e > 3.0, 
the agreement is good although for lower values, there is 
deviation. For lower values of e we need larger values of 
A to obtain a better fit (not shown in the figure). 

By fixing the coupling e = 0.8, we vary the reset time 
r in the range (0,2.0) in steps of .01 and in each case 
the maximum value of t 2 for stability of synchronization 
is calculated using the same criterion. The results are 
shown in Fig. [5] The r 2mQa: obtained from theory and 
shown in Fig. [3] are reproduced here for comparison. The 
numerical values are found to support the results of the 
theoretical analysis very well. We note that here the A 
dependence cancels out and hence the agreement with 
the theory is much better than that for the r 2 — e plots. 

We consider next two Lorenz systems, where the x- 
system is given by 

x\ = a(x 2 - xi) 

X 2 = CXi — X 2 — X1X3 (31) 

X3 = —bx 3 + xix 2 

This is coupled to an identical y-system using the same 
scheme. Choosing parameter values for chaotic Lorenz as 
a = 10.0, b = 8/3 and c = 28.0, the analysis is repeated 
as in the case of Rossler. Here the time step chosen is 
0.001 and r = 0.05. The r 2 values are varied in the range 
(0,0.1) and the stability limits of e isolated. The results 
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FIG. 5: The maximum values of T2 for which anticipatory 
synchronization is stable in two coupled Rossler systems is 
shown as a function of the rest time r. The solid and dotted 
lines are the similar values from theory reproduced from Fig. [3] 
for comparison. 
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I 1 1 1 1 1 

0.02 0.04 0.06 0.08 0.1 

r 2 

FIG. 6: The maximum values of coupling e for two coupled 
Lorenz systems as a function of the anticipatory time ti . The 
solid curve is the limiting curve from theory same as in Fig. [2] 
for a value of A = 0.0. 



V. CONCLUSION 

We introduce a new coupling scheme with varying time 
delay for synchronization of two systems with delay or 
anticipation. The scheme has the advantage that syn- 
chronization can be achieved with intermittent informa- 
tion from the driver in intervals of reset that can be pre- 
fixed. This also makes a detailed stability analysis an- 
alytically possible because the error dynamics becomes 
discrete. By assuming an average effective Lyapunov 
exponent A, the stability regions and limits of stability 
in the parameters of coupling strength and anticipation 
time are worked out for specific cases. We demonstrate 
the method by numerical simulations in two standard 
systems, Rossler and Lorenz. The general features of the 
stability region in parameter space match with the theo- 
retical stability analysis, but more precise matching with 
the numerical data is not possible. This is understand- 
able since in the analytical calculations / is replaced by 
an effective A and also coupling in all components of x 
and y is assumed while in numerical calculations only 
one component is coupled. The agreement between the 
theory and numerical data is reasonably good for the 
T2max — t plot, since the A dependence cancels out. 

The availability of three new time scales in the dynam- 
ics is suggestive of potent applications especially in secure 
communication. We propose that this technique will be 
especially successful with a bichannel transmission [3§| 
where one channel, that is part of the state space of the 
chaotic transmitter (driver), is used to synchronize with 
the receiver (response) and the other forms the message 
along with the chaotic signal from a different part of the 
state space of the driver. Here since the encrypted in- 
formation or cipher text is not used as the synchronizing 
signal, it can be made really complex and secure. In 
this context our method of synchronization has the def- 
inite advantage that the synchronization channel need 
be transmitted only at intervals fixed by the reset time 
which itself forms part of the key space. This leads to 
bandwidth savings and requirement of noise free channel 
for short times at intervals. Moreover, the enhancement 
in the dimensionality of the key space leads to increase 
in security. The stability analysis reported in this paper 
along with the numerical simulations for standard sys- 
tems helps to fix the accessible regions of the key space 
for better key management. This is being worked out 
and will be published elsewhere. 
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APPENDIX A: CASE < r 2 < r 

Here we analyze the eigenvalues /i± of the map (fTTj) 
given by Eq. (fT5|) to obtain the stability conditions. 
These stability conditions are shown in Fig. The syn- 
chronized state, A = 0, is stable if both the eigenvalues 
satisfy \[i± \ < 1. 



1. r 2 = 

For this case 6 = 0. Hence, 2-d map in Eq. (fTTj) be- 
comes a 1-d map given by 

A m+i = fiA m (Al) 

where /i = a (l — (1 — ^)e Ar ). Fig. [7^ shows /i as a 
function oil /a. The fixed point A = is stable provided 
\fx\ < 1. This gives the following limits on a for the 
stability of the synchronized state. 

1 + er XT 



1 < a < 



1-e 



-Xt 



(A2) 



2. < t 2 < r 

For < t 2 < r, the eigenvalues /x± (Eq. {13])) display 
a rich behavior. Three different scenarios are possible. 
These are shown in Figs. [7]b, [7J: and[7Ji which show /i as 
a function of 1/a. To determine the limits of stability of 
the solution A = we consider the following case. 



a. n = 1 (a 2 - 46 > 0, a > 0) 



Putting = 1 in Eq. (fT3j) , we get 2 = a ± Va 2 - 4b. 
This reduces to 



1 = a- b 



(A3) 



Using the expressions (ITUl) for a and b, we get the lower 
limit on the stability as 



ai = 1. 



b. u. = -1 (a 2 - 46 > 0,a < 0) 



(A4) 



Putting fi = -1 in Eq. (|T3j) , we get -2 = a± \/a 2 - 46 
which reduces to 



l+a+6=0 
Using the expressions (|10p . we get 

e Ar + 1 



(A5) 



(A6) 



The above expression gives the upper limit of stability 
for smaller values of r 2 . For larger values of t%, we use 
the condition — 1 which is considered in the next 
subsection. 
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FIG. 7: This figure shows the eigenvalues /i as a func- 
tion of 1/a. (a) Xt — 0.25 and Ar 2 = 0. Here, fi = 
a (1 - (1 - ^)e Ar ) (see Eq. (TSTj). (b) Ar = 0.25 and Ar 2 = 
0.02. The largest /j, (solid line) starts from a value greater 
than one for 1/a > 1, crosses 1 at 1/a = 1, and continues 
till it meets the dashed line from below (Here a 2 — 46 = 0). 
Then fi becomes complex and the dotted line shows the mag- 
nitude \fj,\. This continues till we have a 2 — 46 = again. This 
point is just above the meeting point of solid and dashed lines 
on the negative side. For smaller values of 1/a, /i again be- 
come real (but now negative) and the largest n in magnitude 
jumps to the dashed line bellow. Hence the stability range is 
from a = 1 till the point where fi = — 1. (c) Ar = 0.25 and 
At 2 = 0.1. This figure is similar to (b), but here the dotted 
line (/J, complex) crosses the magnitude one before jumping 
to the negative value. Hence the stability range is now from 
a = 1 till the point where the dotted line crosses one or the 
complex fi has magnitude one. Crossover from the behavior 
(b) to (c) occurs at the peak value, a p , as seen in Figure [2] 
(d) Ar = 1.0 and Ar 2 = 0.7. Here, r 2 > T 2max = 0.458.... 
Hence, the largest \fi\ is always greater than one. 



c. \fj,\ = 1 (a 2 — 46 < 0, /i complex) 

Putting \fi\ = 1 (/i complex) in Eq. (fr3"|) , we get 1 = 

\ \/a 2 — (a 2 ~ 46) which reduces to 

6 = 1 (A7) 

Substituting from Eq. (|10b|b we get the upper limit on 
the stability as 

Ou = — (A8) 



1 



The above expression can also be used to determine the 
maximum T2max for a given r. This happens when there 
is always an eigenvalue with magnitude greater than one, 
i.e. when a; = a u = 1. From Eq. (|A8[) we get the 
following expression. 

-Ar\ 



Ar 2li 



ln(l 



(A9) 
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Note that a u = 1 in Eq. (|A6[) gives the same Ti max as 
in Eq. (|A9|) . Fig. ([7H) shows \i as a function of 1/q for 
7"2 > T- 2maa ; where the synchronized state is not stable. 



d. Peak 

The peak value a p is given by the intersection of 
Eqs. (|A6[) and (|A8|) and leads to the conditions 



6=1, and a + 2 = 0. 
From & = 1, we have 



-At 2 



1 



= 1 e 

a 



-At 



(A10) 



(All) 



Substituting this in a + 2 = 0, we get 

e Ar + 3 



1 



(A12) 



The corresponding t 2 value is given by, 

\t 2p = At + ln(e AT + 3) - 2 ln(e Ar + 1) 



Eliminating r from Eqs. (|A12j) and (|A13|) gives Eq. (TIT)) . 



T2 



This is a simple case where a and 6 in Eqs. ([10]) re- 
duce to a = c = e AT , b = d = a(e Ar — 1) The cases 
Appendix A2a and A2c in the above subsection are ap- 
plicable and hence the stability condition for A = is 



1 < a < 



1 



,At 



(A14) 



The eigenvalue equation is 



H n+1 - CV l t> 



0. 



(B3) 



where n > 1. 

The following general conclusions can be arrived at 
using Gersgorin discs. There is one disc with center at 
c and radius d and n discs with center at and radius 
1. All the eigenvalues lie within these discs. For a < 1, 
c > d+l. Hence, the disc with center at c is disjoint from 
the other discs. Thus one root which lies in this disc, 
must always have magnitude greater than one. Hence, 
the lower limit of stability is ai = 1. 

For n — 1, Eq. (|B3j) becomes a quadratic equation. 
This is discussed in Appendix A3. 



Eq. (|B3[) becomes a cubic equation. 

^i 3 - cy? + d = 

At the upper stability limit, // is complex with 
(A13) Thus, fi 2 — 2 cos(#)/i+ 1 is a factor where n = e ±li 



this condition we get the relation 
d 2 + cd = 1 
This gives a quadratic equation in a. 

(e Ar - l) 2 a 2 + e AT (e AT - l)a - 1 = 0. 
Using the correct solution the stability range is 



1 < a < 



— ( \/e 2XT + 4 - e XT ) 

2(e A --l) V Ve +4 ) 



The maximum T2 m ax is obtained in the next 
tion. 



(B4) 

Ml = 1. 
. Using 



(B5) 
(B6) 

(B7) 
subsec- 



APPENDIX B: CASE r 2 = nr 

This corresponds to the case r 2 = in Section IIIB. 
Using Eq. (f2"Tj) . we get the following recursion relation 
(note that n=k+l gives the correct correspondence) 



A,. 



hn+ i = e AT A m+ „ - a(e XT - l)A m (Bla) 
= cA m+n -dA m (Bib) 



where c 



d At 



and d 



ale 



At 



1). This leads to an 



(ri + 1) dimensional map. This map can also be obtained 
directly from the solution ([6]) noting that for t = mr and 
t = (m + 1)t we get A m+ „ and A m+n+ i respectively. In 
matrix form 



/ A„+i \ 
A„ 
A„_i 

V Ai J 



( c 
' 1 
1 



-d\ 





Vo o ••• i o / 



A n _i 
A„o 



V A / 



(B2) 



2. Any n 

For a general n it is not possible to obtain the stability 
range for the synchronized solution. 

It is easy to see that the maximum Timax is obtained 
if there are two degenerate eigenvalues of Eq. (|B3[) equal 
to one at a = 1. This is possible if c = 1/n, d = 1 + l/ra. 
Using the explicit form of c or d and rt = T2 m ax/T 1 we 
get Eq. (J30I) for T 2ma x- 



APPENDIX C: NEGATIVE A 

If Lyapunov exponent A is negative, then Eq. ([5]) can 
be written as 



A = -|A|A-e£ X ( 



mr, (m+ 1) r) 



(CI) 



m=0 



The analysis is similar to that for positive A. Here, we 
summarize the results. 
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where we define a = e/|A| as the normalized dimension- 2 - 

less coupling constant. 



a. T2 — 

For r 2 = 0, b — 0. The stability limits for the synchro- 
nized state are (see Eq. (IA2|) ) 



0.5 



1 < a < 



1 



,-|A|- 



1 - e" 



Alr 



(C4) 



6. < t 2 < r 

In this case Eq (|C2|) leads to a 2-d map as for the 
positive A. The eigenvalue equation and the solutions 
are same as Eqs. (| 12[) and (fl3|) with a and b defined by 
Eqs. (IC31) . 

The lower stability limit is always ai — — 1. For smaller 
values of r 2 (< T 2p ), the upper limit of stability is given 
by (see Eq. (fTi|) ) 



1 



1 + e -l A l r - 2el A K r 2- r ) 



(C5) 



1 

T2/T 



1.5 



FIG. 8: The stability region of the synchronized state in the 
|A|t2 — a plane. The solid line is for At = —0.25 < and 
the dashed line is for At = —0.5. The lower limit of stability 
is ai — —1 (dotted line). For smaller values of T2 < T2 P , 
the upper limit of stability is given by (Eq. ([C5j) while for 
larger values of T2 (t2 P < T2 < r) it is given by (Eq. ([C6j). 
The peak values are (0.072 . . . /0.25 . . . = 0.29, 17.083 . . .) for 
Ar = -0.25 and (0.165 ... /0. 5 = 0.33 9.16 .. .)for Ar = 
-0.5 (see Eqs. JCZl and JCjl]}). F or r < r 2 < 2r the upper 
limit of stability is given by Eq. (|C11|I . We note that the 
stability limits have a similar behavior to that of Fig. [2] for 
positive A. 



c. T2 = r 



For t 2 = t, the stability range is 
- 1 < a < ' 



1 



Air 



(C9) 



while for larger values of r 2 (r 2p < t 2 < r) it is given by 
(Eq. m - 



2. t 2 > r 



3 |A|r 2 _ I 



(C6) 



It is interesting to note that for very large values of the 
coupling constant the synchronized state is unstable. The 
maximum value of a p is given by the intersection of the 
two curves ([C5]) and (IC6l). 



3e 



At 



1 



(C7) 



The corresponding r 2p value is determined by the relation 
|A|r 2p = |A|T + 21n(l + e-l A l r )-ln(3 + e-l A l r ) (C8) 



Let r 2 = kr + r 2 , k = 0, 1, . . . where r 2 < r as for the 
case of positive A. Eq. (|C1[) leads to a map of dimension 
k + 2. The eigenvalue equation is 



/' 



k+2 „,,k+l 



c// +i + 61/i + 6 = 



(C10) 



where c = e l A l r , b\ = a(l — e^^ 7 " 2 T ' > ) and 6 = 
Q, e -|A|r ^ e |A|r 2 _ ^) p or fc = 0, we recover the case 
< r 2 < r. 

For fe = 1, i.e. r < r 2 < 2r, we have a 3d-map. The 
lower stability limit is a; = — 1 . The upper stability limit 
is 



<<„ = ^-(«i + \l a \ +4a 2 ) 



(Cll) 
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0.2 0.4 0.6 0.8 1 

FIG. 9: The maximum values of coupling e for two coupled 
Rossler systems in the periodic region. The solid curve is for 
the values from theory reproduced from Fig. [8] for a value of 
A = 0.0. The agreement is better with the curves in Fig. [2] 
(dotted line) for a value of A = 0.6. 

where a x = e~ 2 l A l T (l - el ) + e |A| (tj-t) _ 1 and fl2 = 

e -2|A|T^|A|r 2 _ ^2 stability limits arc plotted m 

Fig.m 

We have done numerical analysis for negative A using 
two Rossler systems in the periodic region for c = 2.2. 
The stability limits for synchronization in the Ti — e plane 
in this case is given in Fig. [9] The solid line is the curve 
from theory with A = 0.0 . The behavior of the numerical 
results in general agrees with the theoretical analysis. 
However, exact fit is not obtained for any negative A. 
Surprisingly, the fit is better for positive A with equations 
Eq. |14|) and Eq. (JUJ) (dotted line). The reason for this 
behavior is not clear. 



T~2max 

The condition for obtaining the maximum value Ti max 
is that ai = a u . For negative A, we have ai — — 1 and a u 
always remains positive. Hence, unlike the case of posi- 
tive A, the condition for obtaining Timax is never satisfied 
and synchronized state is possible for any t-i or T2 m ax is 
infinite. 



4. T2 = nr 

This corresponds to the case t 2 = of Appendix C2. 
The eigenvalue equation is (see Eq. (|B3|) ) 

H n+1 -cv n + d = 0, (C12) 
where n > 1 and c = e~' A ' r and d — a(l — e~' A ' r ). 

The following general conclusions can be arrived at 
using the Gersgorin discs. There is one disc with center at 
c and radius \d\ and n discs with center at and radius 1. 
All the eigenvalues lie within these discs. For a < 1, d < 
(1 — c). Since c < 1, the disc with center at c lies within 
the circle = 1. Hence, all the roots of Eq. (|C12I) have 
magnitude less than one and the synchronized state is 
stable. Thus for any n there will be range of a values for 
which the synchronized state is stable. This supports the 
conclusion reached in the previous subsection (Appendix 
C3) that T2max is infinite. 

For n = 1, Eq. (|C12[) becomes a quadratic equation. 
This is discussed in Appendix CI. 

For n = 2 we have a cubic equation. The stability 
range is 

2(l-e-l^)V V J 

(C13) 
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